A LOCAL PARABOLIC MONOTONICITY FORMULA ON 
RIEMANNIAN MANIFOLDS 



EDUARDO V. TEIXEIRA AND LEI ZHANG 

Abstract. In this article we establish a local parabolic almost mono- 
tonicity formula for two phase free boundary problems on Riemannian 
manifolds, which is an extension of a work of Edquist-Petrosyan. 



1. Introduction 

In the theory of two-phase free boundary problems, it is weh established 
that regularity of the interface is closely related to asymptotic behavior of 
solution near the free boundary. In 1984 Alt, Caffarelli and Friedman [2] 
established a monotonicity formula to describe the interaction of the two 
pieces of the solution on each side of the free boundary. This formula has 
been extremely powerful in the regularity theory and it reads as follows: Let 
ui,U2 be two non-negative continuous functions in Bi (the unit ball in M") 
such that Aui > (i = 1, 2) are satisfied in distribution. Suppose ui -^2 = 
and ui{0) = U2{0) = 0, then 

J Br FI JBr \A 

is monotone non-decreasing for < r < 1. 

There have been different extensions of the theorem of Alt- Caff arelli- 
Friedman under different contexts. For example, Caffarelli |3j established 
a monotonicity formula for variable coefficient operators, Friedman-Liu |14j 
have an extension for eigenvalue problems. Another important extension 
has been achieved by Caffarelli-Jerison-Kenig [TO] who replace Auj > by 
Aitj > —1 [i = 1,2). Under this new assumption they prove that (/)(r) is 
uniformly bounded for < r < ^. This is called an "almost monotonicity 
formula" . 

Even though there is no monotonicity in the Caffarelli-Jerison-Kenig for- 
mula, it does provide a control of |Vui(0)|-|Vu2(0)| ifni,U2 are both smooth 
at 0. For many free boundary problems, the control of |Viii(0)| • |Vu2(0)| 
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usually leads to important regularity results. Moreover, for some real life 
problems such as the Prandtl-Batchelor problem [U [31 [5] ) and some classical 
problems ( e.g. see Shahgholian the equations may be inhomogeneous 
and we may not have Auj > {i = 1,2) on each side of the free boundary. 
The "almost monotonicity formula" of Caffarelli-Jerison-Kenig is particu- 
larly useful in these situations and has provided a theoretical basis for the 
regularity theory of many new problems (see for example [10^ [TB] ) . 

For two-phase parabolic free boundary problems, Caffarelli [9] established 
a monotonicity formula for two sub-caloric functions: Let ui,U2 satisfy Auj — 



dm > in R'^ X (-1,0) {i = 1,2), uiU2 = 0, ni(0,0) = ^2(0,0) = 0. Let 



0(^) = ~ / / \Vui\'^Go{x, —s)dsdx / / \Vu2fGo{x, —s)dsdx 



is monotone non-decreasing provided that ui,U2 have reasonable growth at 
infinity. 

Clearly this monotonicity formula for parabolic free boundary problems is 
in correspondence with the Alt-Caffarelli- Friedman formula. Later, Edquist- 
Petrosyan [13j derived the " almost monotonicity formula" for ui , U2 satisfy- 
ing Aui — dfUi > — 1 instead of being sub-caloric. Similar to the Caffarelli- 
Jerison-Kenig formula, Edquist-Petrosyan proved the bound of 0(r) for 
< r < i. 

A common feature of all monotonicity and almost monotonicity formulas 
aforementioned is that they are all designed for problems within Euclidean 
spaces. Prom theoretical and application viewpoints it is natural to con- 
sider some free boundary problems on Riemannian manifolds. Indeed, it 
has been pointed out by Caffarelli and Salsa in [12j that the tools devel- 
oped for free boundary problems on Euclidean spaces should have their 
counterparts for free boundary problems on manifolds (page ix of the intro- 
duction). The analogs of Alt-Caffarelli- Friedman monotonicity formula and 
Caffarelli-Jerison-Kenig almost monotonicity formula have been developed 
for the Laplace-Beltrami operator by the authors in [17j. The purpose of 
this article is to derive the analogue of Edquist-Petrosyan formula on Rie- 
mannian manifolds. In forthcoming works we shall use these formulas to 
discuss free boundary problems on Riemannian manifolds. 

Let {M,g) be a Riemannian manifold of dimension n > 2, let B(p, Sp) be 
a geodesic ball around p with radius 5p = min{l, inj^} (inj^ is the injectivity 
radius at p). We shall use the following cut-off function x supported in 




■ n 



then 



B{p,5p): 



X = 1, in i?(p, X = 0, in B(p,5p/2). 
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Let Rm denote the curvature tensor. Assume 

(1.3) \Rm\ + \VgRm\ < A. 

In this article, if we do not mention the dependence of a given constant, it 
is imphed that this constant is either universal or depends only on n, A and 
X- 

Let Q-{p) = Bp{r) x (-r^O) for < r < 6p/2. If n+,n_ G Hl,{Qs^{p)), 
we define w± as w± = u±x and let 

(^(r) = \ [[ \VgW+\^G{x,-s)dVgds f f \V gW-\^G{x, -s)dVgds 

^ JJSrip) JJSr{p) 

1 d{x,p)'^ r) / 

where G{x,t) = j^^-^e Sr = B{p,6p) x {—r'',0), dVg = y/det{g)dx. 

Our main result is: 

Theorem 1.1. Let u± be non-negative continuous functions in Qj (p) that 
satisfy u+ ■ U- = and 

{Ag - dt)u± > -1, in Qp{l) 
in the weak sense. Then u± € Hl^^{Qj^{p)) and there exists C > such that 

(j){r) < C(l + [[ u\dVgds + [[ u^_dVgdsf 

for all r G (0, \5p). 

The proof of Theorem 11.11 is along the lines of [10] and [13] . However, 
since equations, integrals and kernels are defined on a Riemannian man- 
ifolds, many perturbation terms have to be properly controlled in differ- 
ent situations. For example, we need to derive a "perturbed" version of 
Beckner-Kenig-Pipher inequality in Proposition 12. 3i As the reader will see, 
it is crucial in our analysis to have a quantitative estimate of the perturba- 
tion of the eigenvalues in this key inequality. In accordance to [10] and [13] , 
assuming Holder continuity of solution, we can infer a more precise control 
of the functional (j). This is the content of Theorem 13.11 we present at the 
end of the paper. With our "perturbed" version of Beckner-Kenig-Pipher 
inequality in hand, the proof of Theorem 13.11 becomes very similar to the 
corresponding theorem in jlSj and is therefore omitted. 

2. The proof of the monotonicity formula 

Let U{x,t) be the heat kernel in the neighborhood of p, then for t small 
we have 

U{x,t) = (47rt)-te-^( j;0,(x)f ) 
(see [15] P109) where are smooth functions of x, and 



det(fir)"4 = 1 + 0{d{x 



,2> 
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Since we consider a neighborhood of the origin, G and U are clearly com- 
parable in this neighborhood. As a consequence, G can be replaced by U in 
Theorem 11.11 In the following we will mainly use U in our proof. 
Let 



^t= ff \VgW±\^U{x,-s)dVgds, 



To prove Theorem 1 1.1 1 it is enough to prove the bound of (p{r) for < r < S 
where (5 is a small constant depending on n, A and x- The bound for r > 6 
is obvious. Therefore in the proof we only focus on the estimates of for 
k sufficiently large. 

We shall prove the following two key propositions: 

Proposition 2.1. There exist Co,Ci > such that for k > kQ{n,A,x), if 
> Co, 

where 6k = Ci{^ + + ^~^'')- 



Proposition 2.2. There exists e G (0,1) such that for k > ko{n,A), if 
6± > Co and 4^A+^, > A+ then A^^, < (1 - e)A^. 

Theorem 11.11 follows from Proposition l2.1l and Proposition [2?2] by standard 
argument in [101 117j . Here we note that we shall always assume u± to be 
smooth, as u± can be mollified to such that necessary inequalities for 
can be obtained first. Therefore, the conclusion for u± can be obtained 
by passing to the limit. This part of the argument is standard and is 
omitted. The interested readers may look into [T7] for reference. 

The following estimate is important for the proof of both Proposition 12. II 
and Proposition 12. 2t 

\V gW±\'^U{x, —s)dVgds 



< CMT^ + CMr^ wi{x,-r')U{xy)dVg)-^ 
(2.1) +1 [ wl{x,-r^)Uix,r^)dVg. 

where Cm depends on ffg^{u^ + u'^)dVgds. 

Proof of (12. Ij) : From {Ag—ds)u± > —1, we have, by standard computation 
{Ag - ds)iwl) > -2w±x - 4lVgn±| • \Vgx\w± - 2u±\Agx\w± + 2\VgW±f. 
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From the above we have 

2 // \VgW±\''di^ < II (Ag - ds){wl)du 

+2 II w±diy + 4 II i\VgU±\\Vgx\ + u±\Agx\)w±diy 
= h + h + h- 

Here we use the notation: du = U{x, —s)dVgds. We shah also use dv^ = 
U{x,-s)dVg. 

Using w± > and integration by parts we obtain 



(2.2) Ii< [ wlix,-r^)diy-'"\ 

JRri 

Note that we used 

hm / w'^{x,— e)U{x,e)dVg = -w;±(0,0) < 0. 

To estimate I2 we use the following equation easy to be verified by direct 
computation: 

{Ag - ds){w± -s)> 2VgU±'VgX + u±AgX. 
For S2 < si < we integrate the above to obtain 

/ iw± - si)dv'' < [ iw±- S2)dv'^ 

(2.3) + r ( f {2\VgU±\\Vgx\+u±\Agx\)dAds. 

Js2 \ JR" / 



' S2 

For > si > — > S2 > — 4r^, since the support of Vx or Ax stays 
away from the origin, it is elementary to obtain 

/ w±{x, si)dv^^ < inf / w±{x, s)du^ 

Jr" se[-4r-2,-r-2] Jjjn 

(2.4) +Ci(A, ny + C2(M, N)r^ . 

where M is the norm of u± on , iV is a large number. Therefore I2 
satisfies 

(2.5) /2<2r2 inf [ w±{x, s)du' + C{N, M)r^. 
If we further use Cauchy's inequality we have 

l2<C{M,A)r^+ inf / w^ix , s)diy' . 
The estimate of I3 is similar. 

(2.6) Is < C(iV)r^ / / {\VgU±\^ + \u±\^) < C{N,M)r^ [ [ u^. 

J J 01 J JQ- 



■5p " "'^ Sp 
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()2.ip can be obtained easily from (12. 21) . (12.50 and p.6p . As a consequence, 
the following estimates also hold: 

(2.7) // \ygw±\^dv <CMr^ + Cm inf / wi{x,s)du'. 



(2.8) / / \VgW±\'^du < Cmt^ + ^ I ^±(^' *)^^- 



w±{y,s) = —w±{ry,r'^s) 



Cm 

In the following, we shall always re-scale w± as follows: 

1 

-2 

u± is understood similarly. Correspondingly we let 

9ij{-) = 9ij{r-) = + 0{r'^)\ ■ ^ 
be the re-scaled metric, dv and dv'^ are defined as 

dv = r^-JJiry, -r^s)dVgds, dv" = r"-U{ry, -r'^s)dVg, 
then u± satisfy 

{A-g - ds)u± > -1. 

We use O-t to represent the set where w± is positive. The corresponding set 
for 'w± is ^l±. 

Lemma 2.1. For r < 6p and s = — | 



(2.9) 



^^^l + 0(r2) r ^^^j^s < 2(1 + 0(r2)) / \V-gW±\^du' 



Proof of Lemma I2.lt We use dvQ to represent the Gauss measure in 
Euclidean spaces: 

1 

dvn = —e 2 dx. 

(27r)2 

We perform two transformations on the three integral terms in (j2.9p : 

|i()-|-|ps, / ujj^dv^ , I \V g'iB±^dv^ 

to reduce them to the Euclidean case. First, using y = (t)i{x) (</)i to be 
determined) we have 

|V,.i(x)P=,-(.)^^ 

Here repeated indices imply summation. Since gij{x) is symmetric and 
gij{x) = 6ij + 0(r^[xp) for \x\ < r^^S2 ( 62 small), we can choose ipi so 
that 

1 

dy 



dx 



9ij{x) = 6ij + 0{r \x\ 
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Consequently 



dxi dxj 

and (|2.1(J|) becomes |Vgu)±(x)p = I^^P- Moreover the Jacobian of 

the mapping is of the order 1 + 0(r-^|yp) for \y\ < d^r with 63 small. With 
this (pi we combine the Jacobian with the heat kernel: 

duy := dvJ^^ = r^U{r(f)~^{y), —r'^s)J^^dVgds. 

Using the definition of (pi we now have (recall that s = — i) 

(2.11) r"C/(r</.-i(y), -r^s)J^,dV-g = -i-^e-^(l + 0{r^)\y\'')dy. 

(27r)2 

We use dvy to denote r'^'-U{rx, —r'^s)J^^dVg. With these notations, the inte- 
gral forms in (|2.9p become 

\w±\y3 = 



w±{(piHy))dvl 
[ wl{x)di?'= [ wl{cp^\y))d9l 

(2.12) / \V-gW^\^dv'= [ \Vw±i<P^\y))\^di?'y. 

The purpose of the second transformation is to make di'y as close to the 
Gauss measure on Euclidean spaces as possible. To this end we write dVy as 

1 \y\^ 
d^l = 7:r\E(i ~{l + A{y))dy 

where 

\A{y)\<Cr\l + \y\f, \DA{y)\<Cr\l + \y\), \y\ < 5:,r~\ 
where 5^ is a small number. The second transformation is defined as follows: 

(2.13) y = z + ^{z) 

where ip satisfies z ■ ip{z) = ln(l + A{z)) for |z| > 1. It is easy to obtain from 
the estimate of A that 



(2.14) \^P{z)\ < Cr^\z\, \D^{z)\<Cr^, l< \z\ < d^r^'^ 

where 64 is a small positive number. Then extend the definition of ^p to Bi 
in such a way that both |^| and \Dip\ are of the order O(r^) in Bi. 
Using ()2.13p and (|2.14p we verify by direct computation that 

(2.15) -}—e-'^\y\\l + A{y))J,_^ = ^e-il^l'(l + ©(r^)). 
(27r) 2 dz (2tt) 2 
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Let f±{z) = w±{(l) ^{y{z))), since the Jacobian Jdy = l + 0(r^), the three 

dz 

integral terms in (|2.9p are of the form (see (j2.12p ): 



wl{x)dv^= [ fliz)d,.^o{^ + 0{r')) 

(2.16) [ \W-gW±\du'= [ \Vf±{z)\^du^{l + 0{r^)). 

JR" JR" 

Note that in the last equality, we used ^ = id+0{r'^) where id is the identity 
matrix. For f± we use the Poincare's inequality on Euclidean spaces (see 

my- 

Lemma |2. II follows from the equation above and (j2.16p . □ 

The following two lemmas have analogues in |10[ [T3] and their proofs 
are similar to their counterparts in [13], we include the proofs here for the 
convenience of the readers. 

Lemma 2.2. Let w be w+ or suppose 

Vgwl'^du = ar'^ < oo, / / iVgwl^du > 
nnSr J JnnSr 256 

Then there exists C > and Ci(n, A, M) > such that if a > C, {^10 Sr\ 
Szldi/ > Cir^ > 0, where $7 is the set on which w is positive. 

Proof of Lemma 12. 2t Let w be u)+ or w_. Let Vt be the set on which w 
is positive. Then the assumptions become 

\Vgw\'^dh' = a, / / |Vg?Dpdz/ > 
nnSi JJnnSi 256 

We want to show that if a is large, 1^^ H \ S*! > ci(n). As a result of 

2 4 

(j2.7p we have 

^<[[_ \V-gw\'^di> <C + C inf / w^{x,s)du'. 
256 JJnnSi/4 sg[-| - j^] Jr" 

Therefore for a large 

(2.17) inf / w'^{x,s)di?' > 



From fjgi \Vg'w\'^diy = a, we see that 



\Vgw\'^di?'' < 16a 
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except on a set of line measure no more than jq. So there exists a set E of 
Une measure at least | on [— 3, — such that 

\Vgw\'^du' < 16a, seE. 

For each s & E, either \ijd\du^ > ^, or \w\du= < 5- In this latter case we apply 
Lemma |2. II and (|2.17p to get > c(n). Recall that we always assume r 

to be small. So in either case there exists c(n) > such that \w\du'' > c{n) 
for all s (z E. Therefore Lemma 12.21 is established by scaling. □ 

Lemma 2.3. Let w, Cl he the same as those in Lemma \2.2[ assume /i € (0, 1) 
such that 

2 4 2 4 

Then there exists A(^) E (0, 1) such that 

\Vgw\'^di? < A // \S7gw\'^du. 

Si J J Si 

1 



Proof of Lemma 12. 3t 



1^2 



\nnSi\Si\du = / \F{s)\ds 

where is the measure of the positive set with respect to the dP**. We 

know that < 1 — ^ in a set S C [— (^)^, —(3)^] with the line measure 

greater than or equal to ^l^i \ Si\du- 

2 4 

Using ()2.9p we have, for small r and s E, that 

/ d^i-^sfdv' <C f \V-gw{-,s)\^dv' 

If JJ^ng \Vgw\'^di' < ^, there is nothing to prove. Suppose this is not 
the case, then by using (j2.7p and the largeness of a, we have 

^<l^j,%,,),, v.. [-1,-1]. 

Specifically for s (z E we have 

^< [ w'^{-,s)dv <C f \V-gw{-,s)\^dv\ \/seE. 
This implies that 

a\E\ 



\VgW\^dV > 

JnnSi\Si 
2 5 



Lemma 12.31 follows easily from the above. □ 

The following proposition makes use of the two transformations used in 
the proof of Lemma 12. 1[ 
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Proposition 2.3. Let Cl^^ C M" be the set where W-\-{-, —1) is positive. 
is understood similarly. There exists C > such that 

|Vg^+(,-l)|2dp-i |V,u)_(-,-l)Pdp-i 

+ I - > 1 _ C'r 

!^^W^{;-lfdD'^ + !^,_W.{-,-lfdU-^ - ■ 

Proof of Proposition 12. 3t 

We make the two transformations as used in the proof of Lemma 12.11 
After the transformations, w±{-, —1) become w±, become We still 
have ^\ n llL = and d\. U = i?(0, 6r~^) for some 5 > small. w± are 
supported in Qj^, respectively. Moreover, by the same estimates as in the 
proof of Lemma 12.11 we have 

(2.18) /_ w±{-,-lfdv^^<{l + Cr'^)f w\dv-^ 
and 

(2.19) / \VgW±{-,-l)\^d9-^ y^l-Cr"^) [ \Vw±\^dvQ'^. 



■ \2 



Recall that dun. ^ = — * dx. 

Beckner-Kenig-Pipher inequality (a proof of which can be found in [llj) 
gives 

(2.20) ,° + , > i. 

Therefore Proposition 12 . 31 follows immediately from (I2.18p . ()2.19p and (I2.20p . 

□ 

3. Proof of Proposition 12.11 and Proposition 12.21 

First we observe that Proposition 12.21 is a direct consequence of Lemma 
12.21 and Lemma 12.31 as in \1Q\ I13j. To prove Proposition 12.11 we let 

v±{y, si) = ■^w±{4''y, 42'^si), {y, si) G n±. 

It is easy to see that 0+ and Cl^ are disjoint subsets of -6(0, 5A^) where 5 > 
is small. Let gij be the scaled metric, dD and dD'^ be the new measures. Let 

Also we set 



r"^ JJSr 



A±{r)= [[ \VgV±\^di), B±{r)= I \V gV±\^U{y,-r^)dVg. 

JJSr JM." 
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We want to show that for | < r < 1, if A± = A±{1) are both large, then 
(3.1) 4>'{r) > -C4>{r){-^ + — ^ + 4-2^) 



where C is independent of k. Once (13. ip is estabhshed, the integration of 
(|3.ip gives 



(3.2) <^(i)<^(i)(i + C5,), 5fc = -^ + ^ + 4-2'= 

4 ' ~' 'A^ 



Then by scaUng, (|3.2p is equivalent to Proposition 12. 1[ 

So we are left with the proof of ()3.ip . By scaling, the case for j < r < 1 
can be treated as r = 1. Then the proof is very similar to the standard one: 

,^'(1) = -4i+i_ + 2B+A^ + 2A+B^. 

If B+ > 2A+ or B- > 2A-, 4>'{1)>0. So we only assume B± < AA±. 
Now we apply (|2.ip to v± to get (using dD~^ = JJ{-, l)dVg) 

A±<C + C{f v'idv-^)^ + \ f i^di)-^ 

JR" 2 J]gn 

(3.3) < C + -|=/S + ^B±. 



Note that we can assume A+ and A_ are both positive, because if, say i}+ = 0, 
we obtain from the first line of (13.31) that 



A+<C, 

which is a contradiction to the largeness of Aj^ . 

From (j3.3p we see that if A+ > 2 or A_ > 2, (j3.ip is established easily. 
Therefore we assume \± < 2. In this case we obtain from (j3.3p that 

(3.4) 2A+i+ < C7 + -^V^ + ^+- 

Aj 



There is a similar equation for 2A_^_. Multiplying A^ to p.4p . ^4+ to the 
corresponding equation, we have, by adding these two equations 



(3.5) 2(A+ + A_)i+i_ < C(l + i_^|± + i+^|^) 

+B+A^ + B^A+. 

Proposition 12.31 gives A+ + A_ > 1 - C4-2fc. Using this in (|3.5p we obtain 
(j3.ip . Proposition 12.11 is established. □ 

Theorem 11.11 follows from Proposition 12.11 and Proposition 12.21 as in [10] 
and fia. 
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If further information on the growth of u± is known near the origin, then 
the behavior of (p{r) can be made more precise. This is the observation in 

[iniiis]. 

Theorem 3.1. Let u±,w±,x be the same as in Theorem suppose in 
addition that 

|U±(X,S)| <Ce(|x|2 + |s|)f 

for {x, s) E Qsp and e G (0, 1]. Then 

0(r) < (1 + p')Hp) + Cmp', 0<r<p<6p/4 
where Cm depends on n, A, M = ||u+||2,2(q^^) + \\u-\\i2(^Qg^),x,e- 

Since the proof of Theorem 13.11 is similar to its analogue in [13] , we leave 
the detail to the interested readers. It is worthwhile to point out here 
that the only difference in the proof comes from the correction term in 
the Beckner-Kenig-Pipher inequality, [6j. The readers can easily see that at 
this point the extra term does not cause further difficulties in the argument. 
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